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SUMMARY 

An efficient, decoupled direct method for calculating the first order 
sensitivity coefficients of homogeneous, batch combustion kinetic rate equa- 
tions Is presented. In this method the ordinary differential equations for the 
sensitivity coefficients are solved separately from, but sequentially with, 
those describing the combustion chemistry. The ordinary differential equations 
for the thermochemical variables are solved using an efficient, Implicit method 
( LSODE) that automatically selects the steplength and order for each solution 
step. The solution procedure for the sensitivity coefficients maintains accu- 
racy and stability by using exactly the same steplengths and numerical approxi- 
mations. The method computes sensitivity coefficients with respect to any 
combination of the Initial values of the thermochemical variables and the three 
rate constant parameters for the chemical reactions. The method Is Illustrated 
by application to several simple problems and, where possible, comparisons are 
made with exact solutions and those obtained by other techniques. 


INTRODUCTION 


The modeling of a homogeneous, batch combustion system requires the solu- 
tion of first order ordinary differential equations (ODE's) for thermochemical 
variables such as composition, temperature, and density. The chemistry Is 
represented by a system of NR simultaneous reactions among NS different spe- 
cies. All chemical reactions considered In this work are gas-phase elementary 
reactions. The jth chemical reaction can be written In the general form 

NS 

S ”u x i 

1=1 

where and are the stoichiometric coefficients of species 1 

(with chemical symbol ) In reaction j as a reactant and as a product, 
respectively. 

The time rate of change of species 1 can be written as (refs. 1 and 2) 
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where Is the mole number of species 1 (l.e., mole 1/g mixture), t 

the time, T the temperature, and p the mixture mass density. The net rate 
of formation (f^) of species 1 due to all forward and reverse reactions Is 
given by 
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where the molar forward (Rj) and reverse (R.j) rates per unit volume for 
reaction j are given by 
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The forward rate constant (kj) Is given by the modified Arrhenius 
expression 


kj = AjTOj exp ( Ej/Rl) (6) 

where Aj, nj, and Ej are constants and R Is the universal gas constant. 
Each reaction may be either reversible (bidirectional) or Irreversible (uni- 
directional). For reversible reactions the reverse rate constants (k_j) are 
calculated from kj and the concentration equilibrium constants (Kj) using 
the principle of detailed balancing (ref. 3) 


k -j - Xj/Kj 


(7) 


where Kj Is a function of temperature alone. Two different types of batch 
reaction problems can be Identified: constant and variable density. For a 
variable density problem, the pressure-versus-tlme profile Is given and the 
ODE's for temperature and density take the form (ref. 2) 
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where p Is the absolute pressure and y, A, B, and D are given by 
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In these equations h^ and c p ^ are the molar specific enthalpy and 
constant-pressure molar specific heat, respectively, of species 1, c p the 
mixture specific heat, y the mixture specific heat ratio, M w the mixture 
mean molar mass, and Q the heat loss rate per unit mass of mixture. The 
thermodynamic properties h^ and c p j are functions of temperature 
alone and are computed by using polynomial equations (ref. 2). A can be 
described as a species production function, and B and D are enthalpy 
production and loss functions, respectively. 


For constant density problems the temperature ODE becomes (ref. 2) 

= T £(y - 1)A - B - Y oj (16) 

and the pressure Is obtained from the Ideal gas law 

p = pRT/Mw (17) 

The problem Is to determine the thermochemical variables at the end of a 
prescribed time Interval, given the Initial conditions and the chemical reac- 
tion mechanism. 


The use of classical methods such as the explicit Runge-Kutta and Adams 
methods to solve the ODE'S arising In combustion chemistry results In prohibi- 
tive amounts of computer time. This Is due to the extremely small steplengths 
required by these methods due to the "stiffness" exhibited by the ODE's 
(refs. 1, 4, and 5). The phenomenon of stiffness In chemical kinetic rate 
equations was first recognized by Curtiss and Hlrschfelder (ref. 6) who devel- 
oped a simple backward differentiation method for handling such equations. 

Since then, many approaches have been proposed for stiff ODE's In general 
(refs. 7 to 12), and chemical kinetic rate equations In particular (refs. 1, 

2, and 10 to 20). In several recent publications (refs. 1, 5, 21, and 22) the 
accuracy and efficiency of many techniques for the solution of stiff ODE's 
arising In combustion chemistry have been examined. These studies showed that 
the packaged code LSODE (ref. 23) Is at present the most efficient and accurate 
algorithm for batch combustion chemistry problems. This code has therefore 
been adopted In the present work. 

In addition to solving the ODE's for the thermochemical variables. It Is 
often necessary to know how sensitive the solution Is to the Initial conditions 
and the chemical reaction mechanism parameters. Such a need arises In the 
development of reaction mechanisms from experimental data (ref. 24). The rate 
constants are often not well known and In general, the experimental data are 
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not sufficiently detailed to accurately estimate the rate constant parameters. 

Sensitivity analysis helps to determine the effects of uncertainties In the 
rate constant parameters and the Initial conditions on the solution, and to 
Identify which variables are sensitive to which parameters. The analysis helps 
Identify unimportant elementary reactions which can be discarded, thereby sim- 
plifying the mechanism and hence the problem. At the same time, reactions that 
need additional experimental study can be Identified. Sensitivity analysis ' 

helps also In the understanding of complex reaction mechanisms by showing which 
parts of the mechanism are Important for a given problem. 

i 

Another motivation for sensitivity analysis Is the Increased Interest In 
the numerical simulation of multidimensional chemically reacting flows. The 
coupling of realistic chemical reaction mechanisms with multidimensional flow 
calculations results In prohibitive amounts of computer time. Sensitivity 
analysis can help reduce the reaction mechanism by Identifying those reactions 
that are Important In computing the variables of Interest. 

The simplest method for estimating the effect of uncertainties In any 
parameter Is to run the simulation program with two different values of the 
parameter. Although such a "brute force" method has been used successfully 
(refs. 25 and 26), It can become very expensive when the number of parameters 
Is large. Many methods (refs. 27 to 39) have been proposed for a more elegant 
and less expensive sensitivity analysis. The commonly used methods Include the 
direct method (DM) (ref. 27), the Fourier amplitude sensitivity test (FAST) 
method (ref. 28), the Green's function method (GFM) (refs. 29 and 30) and Its 
variants (refs. 32 and 33), the analytically Integrated Magnus modification of 
the Green's function method (GFM/AIM) (refs. 30 and 34), and the decoupled 
direct method (DDM) (refs. 37 to 39). For stiff problems of the type examined 
In the present work, the DDM has shown greater efficiency and stability, with 
equal or better accuracy than the GFM and the GFM/AIM methods (ref. 37). 

Advantages of the DDM method over other methods are discussed by Dunker 
(ref. 37). 

Dunker (ref. 40) has recently developed a very efficient computer program 
CHEMDDM which couples the DDM method for sensitivity analysis with the code 
LSODE for solution of the model equations. However, the code Is restricted to 
constant density and constant temperature problems. This restriction to con- 
stant temperature problems appears to be a common feature among all sensitivity 
codes developed to date. However, combustion kinetics Is characterized by a 
narrow region of rapidly varying temperature. There Is therefore a need to 
Incorporate sensitivity computations Into a general kinetics code for nonlso- 
thermal problems. In the present work, the DDM method developed by Dunker 
(refs. 37 and 40) Is extended to calculate the first order sensitivity coeffi- 
cients ay-^ /3a, where y^ Is the 1th thermochemical variable and a Is a 
parameter of Interest (either a rate constant parameter or an Initial value), 
for combustion ( 1 . e . , nonl sothermal ) kinetic rate equations. 


DECOUPLED DIRECT METHOD 

The ODE's for batch combustion chemistry presented In the previous section 
can be generalized as 
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where ^ Is the solution vector with NS + 2 components (c^;1=l,NS, T, p) 
and A, n, and E are constant vectors (containing the rate constant param- 
eters, eq. 6), each of which has NR components. 


At present, the most efficient method for solving eq. (18) Is the backward 
difference formula (BDF) Included In the packaged code LSOOE (ref. 1). This 
variable-step, variable-order method computes approximations Y n [= Y(t n )] 
to the exact solution y(t n ) using linear multistep formulas of the type 


* VoLi 


(19) 


where h n (= t n - t n _]) Is the stepslze, q n Is the order of the numerical 

approximation, Y n [= f.(Yjj )] Is the approximation of the exact derivative 

i(t n )[= !(Y(tn)» )]» and the a's and B 0 are coefficients associated with 

the order q n . 

A variety of Iteration techniques Is Included In LSOOE to solve 
equation (19). For combustion kinetics problems the modified Newton Iteration 
procedure Is the most efficient (ref. 1) and Is given by 
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where I Is the Identity matrix, J the Jacobian matrix with element 
•1^. = 3Y /3Y , the superscripts (m) and (m+1) denote the Iteration numbers, 

J «<"» f/ Y (m) 

-n " LV1 n ’ • 


and 


). For each solution step the code automatically 


selects the stepslze and order to minimize the computational work while keeping 
the estimated local error within a user-specified error tolerance. 


The difference equation for the first-order sensitivity coefficient 
n j^= -j^n^ = (aY/a*^)^ j where Is an Initial value or a rate constant 

parameter, can be obta1ned J by differentiating equation (19) with respect to 
nj. The difference equation can also be derived by first differentiating 
equation (18) with respect to nj and then applying the BDF to the resulting 
equation. This gives 


3Y 
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where 3Y n /dnj accounts for any explicit dependence of Y r 
When ni Is an Initial value this term vanishes, and when nj 
constant parameter. It can be obtained from equation (18). Tn( 
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of Sj Is the jth column of the Identity matrix If nj Is the jth element of 
y 0 - If. however, nj Is a rate constant parameter, Sj(t = 0) Is equal to the 
null vector. 

Equations (20) and (21) show the similarity between the model and sensi- 
tivity equations. The DOM method exploits this similarity by alternating the 
solution of equation (19) with that of equation (21). More specifically, the 
solution procedure Is as follows. For any step the solution Y Is advanced 
from t n _] to t n using a standard predictor and the corrector formula 
equation (20). The solution In at the new time step Is then used In equa- 
tion (21) to update the sensitivity coefficients. The latter process does not 
require either a predictor or an Iterative procedure. In solving equation (21) 
for the sensitivity coefficients the method uses exactly the same stepslze and 
order as those used for equation (19). This Implies that the error control In 
the solution of equation (21) Is determined by the error control In the solu- 
tion of equation (19). As discussed In Dunker (ref. 37) the sensitivity coef- 
ficients calculated from equation (21) are the exact sensitivity of Y with 
respect to nj, apart from computer roundoff errors. 

At each solution step, equation (21) must be solved as many times as the 
number of parameters with respect to which the solution sensitivity Is 
required. However, since the matrix (I - h n B 0 J) Is Independent of the sensi- 
tivity solutions, It has to be LU-decomposed only once, irrespective of the 
number of sensitivity parameters. Hence, although the calculation of the sen- 
sitivities with respect to the first parameter may require considerable work 
to form the Iteration matrix, perform Its LU-decomposItlon, and solve 
equation (21), the evaluation of the sensitivity coefficients with respect to 
the second and subsequent parameters Is significantly less expensive. 

To reduce the computational work associated with the calculation and 
decomposition of the Iteration matrix (I - h n B 0 J) In equation (20), this 
matrix Is not updated at every Iteration In LSODE. For additional savings It 
Is updated only when the solution to equation (20) does not converge. Hence 
the Iteration matrix Is only accurate enough for the Iterations to converge and 
the same matrix may be used over many steps. However, to maintain accuracy In 
the computed S j , the matrix In equation (21) must be recomputed at every 
time step unless J changes slowly (refs. 37 and 39). The updating of this 
matrix at every time step obviates the need to Iterate for S-j and specify 
separate error tolerances for the sensitivity calculations. Notice that 
because the same Jacobian matrix J Is required In both equations (20) 
and (21) no additional programming Is required by the DDM method for either the 
calculation of J or the LU-decomposItlon of the matrix (I - h n B 0 J) In 
equation (21). 


ILLUSTRATIVE EXAMPLES 

The DDM method of sensitivity analysis together with a modified version 
of the packaged code LSODE (ref. 23) were Incorporated Into an existing general 
chemical kinetics computer code GCKP84 (ref. 2). The sensitivity subroutines 
were adapted from the code CHEMDDM (ref. 40). The new general kinetics code 
GCKP87 (ref. 41) has been designed to treat a variety of reaction problems 
Including sensitivity analysis of batch combustion kinetics equations. At the 
user's option the code computes the first order sensitivity coefficients 
dyi/3nj where y^ Is any dependent variable, the species mole numbers. 
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temperature, density, and pressure. The n \ ' s are the Input parameter 
values. This Includes the Initial values of the species mole numbers, tempera- 
ture and density, as well as the values of the rate constant parameters A j , 
nj, and Ej for the jth reaction (eq. 6). The code does not compute the 
coefficient ay^/akj used In the constant temperature computations of other 
Investigators, because the rate constant kj Is a function of the temperature 
(eq. 6). The code GCKP87 also computes the sensitivity coefficients of the 
temporal derivatives of the dependent variables ay^/anj for all the y^'s 
and nj's defined above. 

A preliminary version of GCKP87 has been tested extensively on a variety 
of problems to ensure Its accuracy and efficiency (ref. 42). To test the sen- 
sitivity computations many problems have been examined (ref. 41), three of 
which are presented below for Illustrative purposes. All calculations were 
performed In double precision on the NASA Lewis Research Center's IBM 370/3033 
computer. 

Test problem 1, taken from Dunker (ref. 37), describes the pyrolysis of 
ethane at a temperature of 923 K. This constant temperature, constant density 
problem consists of 5 Irreversible reactions among 7 species. The reaction 
mechanism and rate constants (at T = 923 K) are given In table I, together with 
the Initial conditions. Although the mechanism Is quite small this problem Is 
very stiff and other direct methods have produced Inaccurate results (ref. 37). 
Because this Is a constant temperature problem the rate constants kj (eq. 6) 
are time Invariant. Therefore It Is only necessary to compute sensitivity 
coefficients with respect to kj Itself and not with respect to the Individ- 
ual rate constant parameters Aj, n j , and Ej. This problem was selected 
as a test for the modifications made to the sensitivity routines adapted from 
CHEMDDM (ref. 40) to ensure that Dunker's results (ref. 37) could be 
duplicated. 
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Normalized sensitivity coefficients 


kj 

ak j 


aln<^ 

alnkj 


calculated at two different times by the codes GCKP87 and CHEMDDM are given In 
table II, together with the results obtained by Dunker (ref. 37). All sensi- 
tivity coefficients presented In this table are with respect to the rate con- 
stant ki . To generate these coefficients with the code GCKP87, sensitivity 
with respect to the preexponentlal constant (Aj In eq. 6) Is specified. For 
constant temperature problems the normalized sensitivity coefficients with 
respect to Aj are Identically equal to the normalized sensitivity coeffi- 
cients with respect to kj. To enable accuracy and efficiency comparisons, 
both GCKP87 and CHEMDDM were run with the same values for the local error 
tolerance parameters required by LSODE . The values specified for these param- 
eters were the same as those used by Dunker (ref. 37): 10~ & for the local 

relative error tolerance (RTOL) and 10- 8 for the local absolute error toler- 
ance (ATOL) for all variables. 


The agreement between the two codes Is excellent at both 1 and 20 sec, 
thereby confirming the reliability of the code GCKP87 . Comparisons of the com- 
putational work required by the two codes showed GCKP87 to be significantly 
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more efficient than CHEMDDM. For 2 output stations (1 and 20 sec), GCKP87 
required 9b steps with 138 functional and 21 Jacobian matrix evaluations to 
complete the problem, whereas CHEMDDM required 140 steps with 193 functional 
and 28 Jacobian matrix evaluations. To solve for the composition and sensitiv- 
ity coefficients of all species with respect to all Initial values and all rate 
constants, GCKP87 required approximately 0.65 sec CPU time, whereas CHEMDDM 
required approximately 1.4 sec. The given execution times do not Include the 
time required for code Initialization, preprocessing of the thermochemical 
data, and Input and output. When the number of output stations was Increased 
to 5 (10- 3 , lO- 2 , 1, 10, and 20 sec), the difference In computational work 
required by the codes was even more marked. GCKP87 required 87 steps, 116 
functional and 19 Jacobian matrix evaluations, and 0.61 sec CPU time. CHEMDDM, 
however, required 153 steps, 207 functional and 27 Jacobian evaluations, and 
1.5 sec CPU time. This variation of the computational work required by LS0DE 
with the specified value for the first output station has been observed previ- 
ously (ref. 5). It Is caused by the procedure used In LS0DE to calculate the 
first stepslze to be attempted for the problem. 

The second example, taken from Hwang (ref. 32), Is also a constant temper- 
ature problem, but It permits a comparison with the exact solution. It Is a 
simple first-order reversible reaction 


A 



B 


o A (0) -- 1000, o B (0) = 1 


k] = 1000, k_! = 1 

which describes a rapidly changing system. The solution to this problem Is 


> 6 <°)) * (v A (0) - k_,o B <0>)* 

<*B(t) = °a(°) * °b( 0) ~ °A(t) 


Vi ♦ k -i) 


Sensitivity coefficients were generated at the same output stations as 
Hwang (ref. 32) and the results are presented In tables III and IV, together 
with the exact solutions. The tolerances used for this problem were 
RTOL = 10-6 and ATOL = 10~ 8 to be consistent with Hwang (ref. 32), whose 
results were generated with RTOL = 10-&. Table III gives the sensitivity 
coefficients of o A with respect to the Initial values. In table IV the 
sensitivity coefficients with respect to the rate constants are presented, 
along with Hwang's results obtained with a scaled Green's function method 
(SGFM) (ref. 32). These tables show the excellent agreement between the GCKP87 
and exact results. Although the GCKP87 and SGFM results agree well, there are 
some discrepancies In the SGFM results at early times. In particular, the SGFM 
sensitivity coefficients with respect to k_-| at 1.5xl0 -4 and lxl0 - 3 sec 
are noticeably Inaccurate In comparison to the GCKP87 results. GCKP87 was also 
significantly more efficient than the SGFM method. To solve for the composi- 
tion and the sensitivity coefficients of both species with respect to both Ini- 
tial values and both rate constants GCKP87 required approximately 0.3 sec on 
the IBM 370/3033 computer. In contrast, Hwang (ref. 31) states that the SGFM 
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method required approximately 8 sec on a CDC-CYBER- 1 72 computer. GCKP87 Is 
seen to be significantly faster than the SGFM method even after accounting for 
the relative slowness of this CDC computer. Hwang (ref. 32) also attempted 
this problem with a direct method (DM) and gives a CPU time of approximately 
7 sec, which Is significantly longer than that required by GCKP87 . Since no 
DM results are given In (ref. 32), an accuracy comparison with GCKP87 Is 
precluded. 

The last example Is also a simple problem for which an analytical solution 
Is known. It Is, however, a nonconstant temperature problem Involving a first- 
order Irreversible reaction 


k 

A -* B 


k = A in exp( -E/Rl ) 
A = l,n = l,E = 0 


®A(0) = 1, ob(0) - 0, T( 0) = 1000 K 


To solve the problem analytically the following simplifying assumptions were 
made: (1) constant pressure, adiabatic reaction, and (2) constant and equal 
specific heats (c p ) for species A and B. The solution Is 


o A (t) = Co A (0)e' ACt / C - \® A (0)(l - e~ ACt ) 


* B (t) = ® A (0) ♦ ® B (0) - <» A (t) 


where 


T(t) = T(0) ♦ x[o A ( 0) - a A “»] 

x ■ Gc/pIV 0 ’ * V 0 ’] 

C = T ( 0 ) ♦ Xo A (0) 


where Q c Is the heat of combustion which dictates the temperature rise due 
to the reaction. 


The analytical and computed sensitivity coefficients with respect to A 
at various times are given In table V. This solution was obtained with values 
of Q c = 5000 cal/mol and c p = 5 cal/mol K which give a 1000 K temperature 
rise when reactant A Is completely converted to product B. Also presented In 
this table Is the mixture temperature which gives an Indication of the extent 
of the reaction. The agreement between the analytical and the computed results 
Is excellent at all levels of reactedness. In table VI the sensitivity coeffi- 
cients of the time derivative of are presented to Illustrate the capa- 
bility of GCKP87 to compute these quantities. The analytical solutions were 
obtained by differentiating the ODE for a & with respect to the rate con- 
stant parameters. The solution with respect to A Is exact, but those with 
respect to n and E were obtained by using the computed values for 3a^/3n, 
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aT/dn, do/y/a E, and dT/at, because these quantities cannot be obtained analy- 
tically. Ihe normalization procedure used for A cannot be used for n and 
E because these parameters can have zero values. The normalization parameters 
given In table VI (1/ln T and -RT for n and E, respectively) produce a 
1 percent change In k. The normalized sensitivity coefficients with respect 
to n or E can therefore be Interpreted as the percent change In 6 ^ due 
to the change In the parameter n or E that produces a 1 percent change In 
k. This Is analogous to the sensitivity coefficient alna^/alnA which 
represents the percent change In 6/\ due to a 1 percent change In A. For 
a constant temperature problem the three normalized sensitivity coefficients 
are Identically equal to one another. Again, the agreement between the analy- 
tical results and those generated by GCKP8/ Is excellent, Illustrating the 
accuracy of this code for sensitivity analysis of nonlsothermal combustion 
kinetic rate equations. Although this problem Is simple, selected because the 
analytical solution exists, GCKP87 has been used successfully on a realistic 
combustion kinetics problem Involving 110 reactions among 36 species (ref. 43). 


CONCLUSIONS 

An efficient decoupled direct method for calculating the first order sen- 
sitivity coefficients of nonlsothermal combustion kinetic rate equations has 
been developed. Sensitivity coefficients of all thermochemical variables and 
their temporal derivatives with respect to any combination cf Initial values 
of dependent variables and the rate constant parameters of the chemical reac- 
tions can be computed. The method was Illustrated by application to both con- 
stant and varying temperature problems. The computed coefficients agreed well 
with both analytical solutions and those obtained with other codes and solution 
methods. The efficiency of the described method compared very favorably with 
that of other sensitivity analysis techniques. 
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TABLE I. - REACHON MECHANISM FOR 
PROBLEM l a 


Reaction 

Rate 


Constant, 


kjb 

(1 ) c 2 h 6 .ch 3 + ch 3 

1.14 (-2) 

(2) CH 3 * C 2 H 6 -> CH 4 , C 2 H 5 

1.19 (6) 

(3) C 2 H 3 -> C ? H 4 f H 

1.57 (3) 

( 4 ) H f C A - H 2 * C 2 H 5 

9.72 (8) 

1 

(5) H * H -> H 2 

6.99 (13) 


a See reference 37. The Initial concen- 
tration of C2H5 Is 5.951 (-6) mol cm-3; 
all other Initial concentrations are 
zero. 

b Un1ts are mol, cm, s, and the tempera 
ture Is 923 K. Numbers In parentheses 
are powers of 10. 


TABLE II. - SENSITIVITY COEFFICIENTS FOR EXAMPLE 1 


Species 

alncr^/dlnki at 

1 .0 sec 

3ln<^/alnki at 

20.0 sec 


[37] 

GCKP87 

CHEMUDM 

[37] 

GCKP87 

CHEMD0M 

C 2 H 6 

-0.044 

-0.043 

-0.042 

■ 0.820 

-0.800 

-0.789 

CH 3 

1 .000 

1 .000 

1 .000 

1 .000 

1.000 

1 .000 

CH 4 

.976 

.977 

.97 7 

.643 

.650 

.655 

C 2 H 5 

C 2 H 4 

H 2 

.662 

.662 

.662 

-.210 

-.191 

-.181 

.681 

.680 

.679 

.323 

.329 

.332 

.602 

.602 

.601 

.221 

.227 

.230 

H 

.478 

.478 

.479 

.090 

.100 

.106 


TABLE III. - SENSITIVITY COEFFICIENTS WI FH RESPEC1 
TO INITIAL VALUES FOR EXAMPLE 2 


[Numbers In parentheses are powers of 10.] 


t. 

sec 

3 <, a/ 3o a(0 ) 

a<JA/3<JB(0) 

EXACT 

GCKP87 

EXACT 

GCKP87 

1 . 5x1 0“ 4 

0.861 

0.861 

1 .39( - 4) 

1 .39(-4) 

10~ 3 

.368 

.368 

6 . 32( - 4) 

6 . 32( - 4) 

10' 2 

1 .04( 3) 

1 . 04 ( - 3) 

9 . 99( - 4) 

9.99(4) 

1 . 5x1 0“ 2 

9.99(-4) 

9 . 99 ( - 4) 

9.99(- 4) 

9 . 99( - 4) 

2 .0x1 0 -2 

9 .99( - 4) 

9 .99( - 4) 

9 . 99 ( - 4) 

9.99(-4) 
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TABLE IV. - SENSITIVITY COEFFICIENTS WITH RESPECT TO RATE 
CONSTANTS FOR EXAMPLE <* 


[Numbers In parentheses are powers of 10.] 


t, 

sec 

dlnc?/\/3lnki 

dlno A /dlnk_i 

EXACT 

GCKP87 

[32] 

EXACT 

GCKP87 

[32] 

1 .5x1 O' 4 

-0 . 1 50 

-0.150 

-0.148 

1 .20(-5) 

1 . 20( - 5) 

1 . 35( -5) 

10~ 3 

- .999 

-.999 

-1 .03 

7 . 19( -4) 

7 . 1 9( -4) 

8 . 1 0( - 4) 

TO' 2 

-1 .390 

-1.386 

-1.37 

.956 

.956 

.958 

1 .5xl0' 2 

-1 .000 

-1.003 

-1.00 

.999 

.999 

.999 

2 . Oxl 0~ 2 

-.999 

-.999 

-.999 

.999 

.999 

.999 


TABLE V. - SENSI1 i V 1 1 Y COEFF ICIENTS WITH RESPEC1 TO RATE PARAMEIER A FOR 

EXAMPLE 3 

[Numbers In parentheses are powers of 10.] 


t. 

sec 

T. 

K 

dlnog/dl nA 

alnog/alnA 

3lnT/3lnA 

EXACT 

GCKP87 

EXACT 

GCKP87 

EXACT 

GCKP87 

10-6 

1 QQ1 

-1 .001 (-3) 

-1 -001 (-3) 

1.000 

1.000 

■ i 

9.990(-4) 

9 . 990( 4) 

10-5 

1010 

i -1 . 01 Of -21 

-1 . 01 0( -2 1 

1.000 

1 .000 

9 . 900( -3) 

9 . 899( -3) 

l 0-4 

1100 

-.110 

-.110 

.993 

.992 

9 . 003( -2) 

8 . 978( -2) 

10-3 

1761 

-1 .762 

-1.755 

.551 

.550 

.238 

.238 

10“ 2 

2000 

-20.00 

-20.09 

8 . 244 ( -8) 

8 . 1 50( -8) 

4 . 1 22( -8) 

4 .075( -8) 


TABLt VI. - SENSITIVITY COEFFICIENTS OF a A FOR EXAMPLE 3 


t, 

sec 

31 na A /3l nA 

( 3ln<j A /3n)/lnT 

-RT ( dl na/y/3E) 

EXACT 

GCKP87 

EXAC1 

GCKP8 / 

EXAC1 

GCKP8 7 

10-& 

1.000 

mm 

1 .000 

mm 


1 .000 

10-5 

1 .000 

1BKQI 

1 .000 


1 .000 

.996 

10-4 

.980 


.980 

liEJ 1 

.979 

.979 

10-3 

-.523 


-.473 

m 

-.912 

-.927 

10-2 

-19.00 


-18.98 

Bi 

-20.08 

-20.05 
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